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ABSTRACT 

We measured the radial velocity curve of the companion of the neutron star X-ray transient XTE J2123-058. Its 
semi-amplitude (K2) of 298.5 ±6.9 km s" 1 is the highest value that has been measured for any neutron star LMXB. 
The high value for K2 is, in part, due to the high binary inclination of the system but may also indicate a high 
neutron star mass. The mass function (J2) of 0.684 ± 0.047 M Q , along with our constraints on the companion's 
spectral type (K5V-K9V) and previous constraints on the inclination, gives a likely range of neutron star masses 
from 1.2 to 1.8M©. We also derive a source distance of 8.5 ±2.5 kpc, indicating that XTE J2 123-058 is unusually 
far, 5 .0 ± 1 .5 kpc, from the Galactic plane. Our measurement of the systemic radial velocity is -94.5 ±5.5 km s -1 , 
which is significantly different from what would be observed if this object corotates with the disk of the Galaxy. 

Subject headings: accretion, accretion disks — X-ray transients: general — stars: individual (XTE J2123-058) 
— stars: neutron stars — X-rays: stars 



1. INTRODUCTION 



The X-ray transient XTE J2123-058 is a neutron star low 
mass X-ray binary (L MXB) that was discovered in 1998 (L evine, 
Swank & Smith, 1998) . We optically identified the transient 
( Tomsick et al., 1998a[) and determi ned the ~ 6 hour binary 
orbital period (Tomsick et al., 1998b). Type I X-ray bursts in- 



dicate that the system contains a neutron star, and the relatively 
short orbital period indicates that the companion is a late-type 
star on or close to the main sequence. A pair of high frequency 
quasi-periodic oscillations was detected, which may indicate 
a rapidly r otating neutron star (|Homan et al., 1999|; Tomsick 
et al., 1999). XTE J2123-058 distinguishes itself from other 
LMXBs by having high Galactic latitude (b = -36°) and also 
a high binary inclination (i ~ 73°, Zurita et al. 2000). Both of 
these properties are advantageous for determining the mass of 
the neutron star from optical observations in quiescence, which 
is the primary goal of this work. 

Such mass measurements are important for constraining neu- 
tron star equations of state and for understanding the evolu- 
tion of neutron star systems. Precise mass measurements have 
been made for millisecond radio pulsars (MSPs, Thorsett & 
Chakrabarty 1999), but these measurements are lacking for 
LMXBs. Since it is theoretically possible to spin-up neutron 
stars in LMXBs to millisecond periods via accretion, a link 
between L MXBs and MSPs has long been suspected (A lpar 
et al., 1982). Although there is substantial evidence to support 
this picture, the prediction that rapidly rotating neutron stars 
in LMXBs should be more massive by 0.1 to 0.5 M Q (Bhat- 
tacharya 1995) than those that have not been spun-up has not 
been tested. 

2. OBSERVATIONS AND DATA REDUCTION 

In 2000 August, we performed spectroscopy of XTE J2123- 
058 with the Echelle Spectrograph and Imager (ESI, Sheinis et 



al. 2000) at Keck Observatory. As shown in Table [l , we took 
eight exposures over two nights covering the 6 hour binary or- 
bit. The cross-dispersed spectra were taken in Echelle mode 
with ten spectral orders falling on a 2048-by-4096 pixel CCD. 
We performed on-chip rebinning to 2048-by-2048 pixels to de- 
crease read noise. Although this instrument is sensitive from 
3900-1 1000A, we only used the data in the 4700-6820A band 
where the dispersion varies from 0.36 to 0.52A per rebinned 
pixel. This is sufficient to over-sample the resolution of the 
spectrograph as we used a 1" slit, providing a spectral resolu- 
tion of 1.0 to 1.5 A FWHM. 

We obtained R and V-band photometry of XTE J2123-058 
at MDM Observatory 10 days before the Keck observations, 
and found values for R and V consistent wi th the quiescent 
magnitudes reported in Tomsick et al. (1999). We measured 
V = 22.68 ±0.11 (corrected for atmospheric extinction but not 
reddening), and we use this value below to determine the source 
distance. Due to the faintness of the source, relatively long ex- 
posure times between 2340 and 2700 seconds were necessary 
for our spectroscopic observations. Even these exposures are a 
relatively small fraction of the orbital period. We obtained M, K 
and G dwarf comparison star spectra for the cross-correlations 
described below. 

We extracted the spectra using the software package "ma- 
kee" which was originally developed for the HIRES instrument 
at Keck, but has been adapted to extract ESI spectra. Using 
makee, we produced optimally extracted spectra for each of the 
ten spectral orders. We obtained CuAr line lamp exposures to 
determine the wavelength solution, and we checked this solu- 
tion for each exposure using OI and Nal night sky lines. Av- 
eraging over all exposures and night sky lines, the difference 
between measured and actual wavelengths is 0.018A. The cal- 
ibration is also stable from exposure-to-exposure with a stan- 
dard deviation of 0.027 A or less for all night sky lines, cor- 
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Table 1 

XTE J 2 123-05 8 Keck Observations 





mjdH 


Exposure Time 


FWHM Seeing 


Exposure 


(days) 


(seconds) 


(arcseconds) 


1 


51759.39442 


2700 


0.87 


2 


51759.42822 


2510 


0.91 


3 


51759.45771 


2510 


1.23 


4 


51759.48768 


2510 


0.79 


5 


51759.51748 


2400 


0.85 


6 


51759.54536 


2340 


0.98 


7 


51759.57288 


2340 


0.86 


8 


51760.35770 


2400 


0.67 



200 



"Modified Julian Date (JD-2400000.5) at the midpoint of the exposure. 



responding to a radial velocity measurement error of less than 
1 .5 km s -1 . The data reduction procedure includes heliocentric 
velocity and atmospheric extinction corrections. We used IRAF 
routines to flux calibrate the spectra and obtai ned exposure s of 
the spectrophotometric standard Hiltner 102 (Stone, 1977) on 
both nights for this purpose. We combined the spectral orders 
after flux calibration and rebinned the spectra to a logarithmic 
wavelength scale with 22.8 km s" 1 pixels. Finally, we dered- 
dened the XTE J21 23-058 spectra using Ay = 0.37 magnitudes 
(|Hynes et al„ 200 1|). 



The quality of the XTE J2123-058 spectra varies signifi- 
cantly between exposures as indicated by the signal-to-noise 
ratios (SNR) given in Table || The mean SNR per pixel in the 
5650-5850Aband is 1.66 in the worst case and 6.63 in the best 
case. One reason for this is that the source shows large orbital 
flux variations as can be seen from the fluxes given in Table |[ 
The flux variations show two peaks per orbit (exposures 5 and 
8), which is expected if the variability is due to ellipsoidal mod- 
ulations. XTE J2123-058 light curves from optical photome- 
try show th at ellipsoidal modulations are indeed present (Z urita 
et al., 2000). The SNR variations are also caused by variable 
seeing as shown in Table [l]. The FWHM seeing for the expo- 
sure taken on the second night is 0".67 compared to an average 
of 0".93 for the exposures on the first night. It should be noted 
that the Table || fluxes are not corrected for seeing. 

3. RESULTS 

3.1. Cross-Correlations 

Our analysis of the spectra focuses on the region from 4700A to 
6820A since the SNR becomes very poor below 4700A and 
night sky lines dominate the spectrum above 6820A. In this 
band, the spectra show features from the accretion disk (Ha 
and H/3 emission lines) and from the companion star. An ex- 
amination of the spectra clearly shows the presence of a wide 
depression centered between 5100 and 5200A. This is the sig- 
nature of a K-type star and is due to the MgH feature at 5180A, 
the TiO band at 4954-5200A and a continuum discontinuity as- 
sociated with the Mgl triplet at 5168-5185A (McClintock & 
Remillard 1990 and references therein). An absorption feature 
is seen near 5900A due to the Nal doublet at 5891-5898A. The 
strength of this feature indicates a late-K spectral type for the 
companion. 

We used the IRAF routine "f xcor" (iTonry & Davis, 19791; F il- 
ippenko, Matheson & Ho, 1995) to cross-correlate the XTE J2123- 
058 spectra with the comparison stars. After masking night 
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FIG. 1. — Radial velocity curve fitted using a sinusoid with the parameters 
given in Table W and folded on the orbital period. Heliocentric velocities are 
shown (with errors increased by 60% over the values given in Tablebn, and the 
dashed line marks the systemic velocity. The data are plotted twice Tor clarity. 



sky lines at 5205, 5579, 5891, 5898 and 6302A and the Ha 
and H/3 emission lines from the source, we used the wave- 
length bands 4890-5190, 5210-5560, 5630-5860, 5985-6135 
and 6640-6820A for the cross-correlations. We also applied 
a bandpass filter from 30 to 800 cycles per spectrum in or- 
der to remove contributions due to spurious variations in the 
continuum and high frequency noise. We calculated cross- 
correlations for GOV, G5V, G9V, K0V, K2V, K4V, K7V, MOV 
and M1V with the exposure 8 spectrum. Highly significant 
peaks with Tonry & Davis (1979) R- values between 8.9 and 
16.3 at roughly the same velocity occur, with the M-type and 
GOV stars giving the weakest correlations and the K7V star 
(HD 88230) giving the strongest correlation. For this reason 
and because of the strong Nal line, we used the K7V star for 
the radial velocity study. 

We obtained unique and significant cross-correlation peaks 
at reasonable velocities (-1000 to 1000 km s" 1 ) for all eight 



Table 2 

Fluxes, Signal-to-Noise Ratios and Cross-Correlation 
Results 



Exposure 


Fluxf' 




snr|' 




R : 


Measured RV 






Fitted RV 






1 


1.24 


1.66 


3.26 


55.3 ±17.6 




63.2 


2 


1.72 


2.66 


4.29 


-191.0±22.7 


-182.4 


3 


1.53 


2.15 


4.54 


-298.1 ±24.6 


-352.8 


4 


2.90 


4.40 


8.88 


-390.1 ±12.2 


-385.0 


5 


3.06 


4.18 


8.82 


-265.2 ± 12.5 


-259.0 


6 


2.48 


2.87 


5.14 


-69.6±17.3 




-58.1 


7 


2.32 


2.41 


3.35 


192.2 ±25.0 


123.6 


8 


5.53 


6.63 


16.32 


191.0±7.4 




193.5 



"Mean flux in the 5650 to 585()A band in units of 10~ 18 erg cnT 2 S" 1 A -1 . 
;, Mean signal-to-noise ratio per pixel in the 5650 to 5850A band. 



Tonry & Davis (1979) R- value 



with 



Tonry & Davis (19791 



Measured heliocentric radial velocity in km s 
errors. 

''Heliocentric radial velocity in km s _1 based on a sinusoidal fit to the data. 



3 




5000 



5500 6000 
Wavelength (angstroms) 



6500 



FIG. 2. — The top spectrum shows the average, Doppler corrected spectrum for XTE J2 123-058 in units of 10 erg cm - 1 s" 1 A' 1 . Ha and H/3 emission lines are 
present at 6564A and 4863A. The spectrum for the K7V comparison star (HD 88230) is shown below and the units are arbitrary. Several matching features between 
XTE J2123-058 and HD 88230 are present including the continuum break near 5200A and the Nal doublet at 5891-5898A. 



exposures. The measured heliocentric velocities range from 
-390 km s" 1 to 192 km s" 1 and are given in Table || The R- 
values vary significantly from 3.3 to 16.3 due to SNR differ- 
ences. Figure p] shows the radial velocity curve. First, we per- 
formed a least-squares fit to the data with a function consisting 
of a constant plus a sinusoid with four free parameters: the sys- 
temic velocity (7), the epoch of inferior conjunction (To), the ra- 
dial velocity semi-amplitude (^2) and the orbital period (P or b)- 
This does not provide a formally acceptable fit (x 2 /v = 12/4) 
so we increased the errors on the individual velocity measure- 
ments by 60% to give a reduced-^ 2 of 1.0. The fit gives a large 
range of possible values for P or \, with a 90% confidence error 
region from 21125 to 22575 seconds. This range contains the 
photometri c orbital period of 21445.5 ±2.3 seconds measured 



previously ([Tomsick et al., 199S ; Uovaisky & Chevalier, 1998 
Zurita et al., 2000|), and we refitted the radial velocity curve af- 



ter fixing P or t, to this value. We expect that the photometric and 
spectroscopic orbital periods are the same since the photometric 
value is based o n partial eclipses seen in the optical light curve 
during outburst (Tomsick et al., 1999; Zurita et al., 2000). 

Table |3| shows the fit parameters with P or b fixed, and this 
fit is shown in Figure [j]. Its semi-amplitude (7^) of 298.5 ± 
6.9 km s" 1 is the highest value that has been measured for any 
neutron star LMXB. The value of -94.5 ±5.5 km s" 1 for 7 is 
consistent with the tentative value of -90 km s" 1 derived from 



emission lines (Hynes et al., 2001). We note that previous mea 



surements of the photometric epoch are too uncertain to make 
a meaningful comparison to our spectroscopic epoch. 

3.2. Average Spectrum and Source Distance 

We produced an average, Doppler corrected spectrum using 
the eight XTE J2123-058 exposures. We first removed pix- 
els contaminated by night sky lines. Next, we multiplied each 
spectrum by a constant based on the fluxes given in Table || to 
establish a common normalization. This is important because 
the wavelength bins near night sky lines do not have contri- 
butions from all eight exposures. We then Doppler corrected 
the spectra using the fitted radial velocity values given in Ta- 
ble ^. Finally, we interpolated the spectra to a common set of 
wavelength bins and calculated the weighted average for each 
wavelength bin. 



The average spectrum is shown in Figure ^| along with the 
spectrum of HD 88230. Gaussian fits to the XTE J2123- 
058 emission lines give equivalent widths (EW) of 15.8 and 
17.6A and FWHM widths of 1443 ± 64 and 1729 ± 133 km s" 1 
for Ha and H/3, respectively. The FWHM widths are signif- 
icantly less for the individual exposures. The continuum fea- 
tures near 5000-5300A are similar to HD 88230, and narrow 
Mgl absorption lines are seen in both spectra near the break in 
the continuum. Matching absorption lines are also present in 
several other parts of the spectrum, including the Nal doublet 
at5891-5898A. 

We fitted the region of the spectrum near the Nal doublet with 
a constant plus two inverted Gaussians for the XTE J2123-058 
average spectrum and also for the comparison stars in order to 
determine the total EW of the Nal doublet. For XTE J2123- 
058, we measure an EW of 3.1 ± 0.6A (after correcting for in- 
terstellar absorption using results from Hynes et al. 2001). For 
the K0V, K2V, K4V, K7V and Ml V comparison stars, we mea- 
sure EW of 1.1, 2.5, 3.6, 6.9 and 7.2A, respectively. Since the 
Ha and H/3 emission lines indicate a significant disk contri- 
bution, these EW values indicate that the spectral type of the 
companion is later than K4V. However, based on the strengths 
of the cross-correlation peaks described above, the companion 
is a K rather than an M-type star. The range of possible spectral 
types is K5V-K9V. 

We det ermined the source dista nce using the method de- 
scribed by McClintock et al. (2001 ). We assume a spectral type 



Table 3 
Spectroscopic Parameters 



Parame[er 


Result 






Epochfj 

Systemic Velocity 
Velocity Semi-Amplitude 
Orbital Period^ 


T = HJD 2451760.0462 ± 0.0013 days 
7 = -94.5 ±5.5 km s~' 
K 2 =298.5 ±6.9 km s" 1 
P ,i, = 21445.5 seconds 



"68% confidence errors are given. 

''Heliocentric Julian Date of the epoch of inferior conjunction. 

c Orbital period dete: 
vaisky & Chevalier, 199 



^ nerl from photometry (Iomsirk et al , 199^; J lo- 
unta et al., 
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FIG. 3. — The neutron star mass (M\) vs. the companion mass (M2). The 
solid line comes from the values of fa and i given in the text, and the dotted 
lines are based on the extremes of the error regions for these parameters. The 
vertical dashed lines mark the likely range for Mi- 



di K7V in the calculation and discuss possible associated un- 
certainties below. The strength of the Nal doublet indicates 
that the fraction of the light coming from the companion at 
5900A is 45%± 8%. This, along with the MDM measurement 
(V = 22.68), gives an effective V-magnitude for the compan- 
ion of 23.56. We calculated the absolute V- magnitude (My ) 
of the companion using the expression given in Popper (1980), 



which depends on the spectral type and the radius of the com- 
panion. Since the companion fills its Roche lobe, its average 
density depends only on the orbital period (F aulkner, Flannery 
& Warner, 1972). From the density and assuming that the com- 
panion mass is 0.4 M w , which is likely correct to within a factor 
of two (Barret, McClintock & Grindlay, 1996), we estimate a 



radius of 0.56 Rq for the companion. This radius and the K7V 
spectral type gives My = 8.42. After accounting for extinction, 
we obtained a source distance of 9.0 kpc. Recalculating with 
different spectral types indicates that the error associated with 
the K5V-K9V range is ±1 kpc. Also, the factor of two uncer- 
tainty in the co mpanion's mass translates to a 25% uncertainty 
in the distance ( |McClintock et al., 2001 ). Finally, there is a 
systematic error due to orbital variations in the V-magnitude. 
The MDM observations were made close to photometric min- 
imum so that the average V-magnitude is somewhat less than 
22.68 . Based on the quiescent light curve shown in Z urita et al. 
(2000), we estimate that the average V-magnitude is close to 
22.55, which decreases the distance by 0.5 kpc. We conclude 
that the distance is 8.5 ±2.5 kpc, indicating that XTE J2123- 
058 is 5.0 ± 1 .5 kpc from the Galactic plane. 

4. DISCUSSION AND CONCLUSIONS 

From K2 and P or b, we derive a mass function of fi = 0.684 ± 
0.047M Q , representing a lower limit on the neutron star mass 
(Mi). The binary inclination (i) and the mass ratio (q = M\ /M2) 
are also necessary to obtain a precise measurement of Mi . For 
XTE J2123-058, i = 73° ±4° wa s determined via m odeling of 
the outburst optical light curves ( Zurita et al., 200C ). Figure^ 
shows Mi vs. M2 using the measured values of fa and i. From 
our constraint on the spectral type of the companion (K5V- 
K9V), it is very likely that M 2 lies between 0.4 and 0.7 M . 
For this range of companion masses, the range of neutron star 
masses is 1.2-1.8Mq. Accurate determination of q, which may 
be possible via measurements of the rotational broadening of 
the secondary or light curve modeling, is critical to improve the 
mass constraint. Finding that the neutron star mass is near the 



upper end of the 1.2-1.8Mq mass range would have important 
implications for neutron star equations of state, and, in general, 
determining the mass has important implications for the evolu- 
tion of neutron star systems. 

The s ource distance we obtain is consistent with previous es- 
timates (Homan et al., 1999; Tomsick et al., 1999), confirming 



that XTE J2123-058 is unusually far from the Galactic plane 
and implying a bolometric luminosity between 3 x 10 37 and 
10 38 erg s" 1 fo r the brightest X-ray b ursts detected during the 
1998 outburst (Tomsick et al., 1999). Given the distance and 



longitude (I = 46.5°) of XTE J2123-058, the systemic radial 
velocity of -94.5 km s" 1 is significantly different from what 
would be observed for an object corotating with the disk of the 
Galaxy. This fact alone does not allow us to distinguish be- 
tween a scenario where the source was ejected from the Galac- 
tic plane and the possibility that XTE J2123-058 has a halo 
origin. 
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